CERN PH TH /2009 067 
MAN/HEP/2009/19 
HERWIG/09/04 
KA-TP-05-2009 
TPJU-1/2009 
MCnet/09/10 



^ I Non-perturbative effects in ttie transverse momentum 

O ■ distribution of electroweak bosons at tfie LHC* 

(N 

Andrzej Siodmok 

Marian Smoluchowski Institute of Physics, Jagiellonian University, 
ul. Reymonta 4, 30-059 Cracow, Poland 

psj I LPNHE, Pierre et Marie Curie Universites Paris VI et Paris VII, 

Tour 33, RdC, 4, pi. Jussieu, 75005 Paris, France 



AND 

Stefan Gieseke 



Oh: 

<~i , Institut fiir Theoretische Physik, 

Universitat Karlsruhe, 76128 Karlsruhe, Germany. 

^ , AND 

■ Michael H. Seymour 

in 

^ . PH Department, CERN, 1211 Geneva 23, Switzerland 

1^ I and School of Physics & Astronomy, University of Manchester, U.K. 

o 
a^ 
o 

The transverse momentum of electroweak bosons in a Drell-Yan process 
is an important quantity for the experimental program at the LHC. The 
, new model of non-perturbative gluon emission in an initial state parton 

^ • shower presented in this note gives a good description of this quantity for 

5^ \ the data taken in previous experiments over a wide range of CM energy. 

The model's prediction for the transverse momentum distribution of Z 
bosons for the LHC is presented and used for a comparison with other 
approaches. 

PACS numbers: 12.39-x, 14.70.Hp, 14.70.Dj, 14.70.Fm, 12.38.Lg 

1. Introduction 

The Drell-Yan process has been widely studied in many past [T] and 
present [21 ^S] hadron collider experiments and played a significant role in 
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the development of our understanding of QCD and electroweak (EW) inter- 
actions, both from the experimental and theoretical point of view. Certainly 
this will also be the case for the LHC experiments, especially because it will 
soon become the unique W and Z-boson production factory which is ex- 
pected to collect 300 million W and 20 million Z events per year of its 
operation at energies ^/s = 14 TeV and the luminosity of 10^'^ cni^'^s^^. 
Among various distributions of W and Z observables, the transverse mo- 
mentum spectrum of vector bosons in a Drell-Yan process is a very useful 
and important quantity for the experimental program at the LHC. In the 
case of W production, the uncertainty in the shape of the spectrum directly 
affects the measurement of the W mass [4J and its mass charge asymmetry 
MiY+ — MiY- [5j. It also helps to understand the signature for Higgs boson 
production at either Tevatron or LHC [6]. Although the experiments mea- 
sure the Z transverse momentum distribution and use this to infer that of 
the W boson, the extent to which the effects are non- universal limits the 
ultimate accuracy of the measurement, unless elaborate tricks, as proposed 
in Ref. [7] are used. For these reasons, it is of utmost importance to predict 
the W and Z observables with as high as possible theoretical precision. The 
sources of uncertainty in the theoretical predictions of observables, such as 
the transverse momentum of electroweak bosons discussed here, are of per- 
turbative and non-perturbative origin. In this short note we will concentate 
on the modelling of the latter in the framework of a backward evolution 
parton shower approach [8] which is widely used in general purpose Monte 
Carlo Generators such as Herwig [9j, Pythia [10] or Sherpa [11]. During the 
parton shower evolution, which terminates at some scale of typical hadron 
mass, the recoil from the emitted gluon^ builds up a transverse momentum 
for the W/Z. In order to fit existing data, the conventional backward evo- 
lution parton shower approach needs to be supplemented by the so-called 
intrinsic (or 'primordial') transverse momentum kj- distribution of partons 
initiating the shower. The physical motivation behind this additional non- 
perturbative ingredient is the Fermi motion of partons within a hadron. 
Therefore, its average value per parton can be estimated based solely on the 
proton size and uncertainty principle to be of the order of 0.3 - 0.5 GeV. 
But the values extracted from data first of all are too large and secondly 
grow with collision energy which cannot be explained by Fermi motion. For 
example, in Herwig+-|- its value grows from = 0.9 GeV, which is needed 
to describe the data taken at the energy ^/s = 62 GeV (experiment R209), 
to 2.1 GeV which, is needed at the Tevatron energies {^/s = 1800 GeV). 
This motivated us to propose a model for backward evolution in which an 
additional non-perturbative component at low transverse momentum pro- 



^ Together with other backward-evolution steps, such as an incoming sea-quark being 
evolved back to an incoming gluon by emitting a corresponding antiquark. 
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vides additional smearing at each step of the evolution. By construction we 
expect more non-perturbative smearing for longer parton shower evolution 
ladder which might cure the problem of dependence on centre of mass en- 
ergy as well as on the size of needed intrinsic smearing, which in our studies 
is kept, according to the Fermi motion argumentation fixed at 0.4 GeV. 

In the following sections we will first briefly describe the model, then 
present how it fits the existing data sets, and at the end of this note we will 
demonstrate the model's predictions for the LHC energies which we use for 
a comparison with other approaches. 



2. Model 

The implementation of transverse momentum production in which non- 
perturbative smearing takes place throughout the perturbative evolution, 
was achieved by a simple modification to an initial-state parton shower 
algorithm. The model was implemented in the framework of Herwig++ 
[12j in which the Sudakov form factor for backward evolution from some 
scale gmax down to q takes the form 



(1) 

with x' = x/z^ for further details cf. Ref. |13j . 

The argument of the strong coupling as in Eq. ([1]) is the transverse 
momentum p_\_ of an emissioro. The cut-off scale represented by p_\_^ is 
needed to avoid divergence of the strong coupling. Below the cut-off scale 
Us is equal to zero and consequently the derivative of the Sudakov form 
factor is equal to zero which translates to zero probability of the gluon 
emission below p^^. Therefore, the two arguments of the Sudakov form 
factor, p^^^^ and p±q, are not the evolution variables but only explicitly 
specify the available phase-space of an emission. 

In order to populate the phase-space below by additional non- 
perturbative emissions we introduce the additional Sudakov form factor 
IS.NP such that 



^ In Herwig++, the argument of Qs is a slightly simplified expression, equal to the 
transverse momentum to the required accuracy, but not exactly. We have tested the 
implementation of our model with this simplified expression and the exact expression 
for transverse momentum, and find very similar results. We therefore use the default 



expression. 
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We achieve this by extending as{p±) into the non-perturbative region using 
the foUowing model 

, . jv{P±), P±<P±o^ /ON 

as[P±) = < (pert). . . (3) 

[a)r '{p±), p±>p±o- 

In order to explore the possibility of a reasonable description of experi- 
mental data, we have studied in a greater detail two simple choices of the 
non-perturbative function (p{p±): flat continuation of as{p_L < P±o), with 
a constant value ipo, as{p_L < p±o) = v'o and a quadratic interpolation 
between the two values asip±o) V'o = fi^)'- 

as{p± < P±o) = ^0 + {as{p±o) - 'fo)p\/p\o ■ 
In both cases our model is determined by two free parameters p±^ and 



3. Data sets and fitting results 

In this section we present some new results of the model which were 
obtained after important improvements of Herwig++'s parton shower, re- 
leased with version 2.3.1 of the program. The main change in the program 
was a fix for a wrongly applied PDF veto in the parton shower q ^ q g 
splittings which, by construction of our model could have influence on pre- 
viously presented results [Hj. Therefore, we have repeated the procedure 
described in detail in [H] and have fitted the two parameters of our model 
to the Drell-Yan data from three experiments: the fixed target p-Cu Fer- 
milab E605 [T5] \/i = 38.8 GeV, CERN ISR p~p colhsions experiment R209 
[16j at ^/s = 62 GeV and CDF Tevatron Run I experiment with energies at 
^/s = 1800 GeV [2]. These experiments cover the whole spectrum of centre 
mass energy for the Drell-Yan process data sets which are interesting for 
our studieaj. 



3.1. Parton-level study 

In the case of purely parton-level shower, with all the light-quark and 
gluon effective masses and cutoffs set to zerqj, with our model for the low- 
scale Us as the only non-perturbative input the fitting procedure gave the 
optimal value for the quadratic extrapolation with as(0) = 0.0 and = 

^ There are more data available but all at even lower CM energies. 

* For technical reasons, it is not possible to set them exactly to zero. However, we have 

confirmed that if they are small enough their precise values become irrelevant and 

have very little effect on the results. 
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Fig. 1. The comparisons of the parton-level results from the non~perturbative 
model with the data from E605 with /bin = 0.88 (left), R209 /bin = 0.76 
(middle) and CDF x^ /bin = 1.0 (right). The Monte Carlo resuhs come from our 
parameter set with ip^ — 0.0, = 0.70 GeV. 



0.7. The resulting \ow-px distributions for the new values are presented on 
top of the data sets in Fig. [TJ The /bin values are a little higher than 
before the parton shower improvements, nevertheless the agreement with 
data remains at a high-level; for all the experimental data sets are below 
or equal one. If we are only interested in the W/Z transverse momentum 
distribution, it is enough to use a parton-level study, however, if one needs 
to simulate fully exclusive events then a hadronization model has to be used. 



3.2. Hadron-level results 



The hadronization model used in Herwig++ requires termination of the 
shower using non-perturbative effective parton masses tuned to e^e~ data. 
Therefore, we performed the same analysis as above but this time with 
restored tuned effective parton masses. In this case the best and most 
stable situation was found for as(0) = 4 and p_Lq = 2.5 GeV, giving the 
per degree of freedom of 0.80 for CDF, 0.66 for R209 and the worst 8.6 for 
E605. We should stress that the used parameter set may not be the optimal 
choice for each experiment or CM energy but rather the best compromise 
between the three experiments. As the fixed target data in our analysis do 
not even include the systematic errors quoted to be around 5-10%, we have 
deliberately put a bit more emphasis on the Tevatron results. 
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Fig. 2. Comparison of the hadron level results from the non-perturbative model 
with data from E605 with x'^/bin = 8.6 (left), R209 x'^/bin = 0.66 (middle) and 
CDF x^ /hin — 0.80 (right). The Monte Carlo results are from our parameter set 
with = 4.0,p_Lo = 2.5 GeV. 



3.3. Remarks 

The first remark is that the new parameter choices for both the parton 
and hadron-level models, are not very different from the ones obtained using 
the old version of Herwig++. Before parton shower improvements our best 
choices were, for the parton-level mode: q;s(0) = 0.0, p^^ = 0.75 GeV, and 
for the hadron-level case: 0^(0) = 3 and p^^ = 3.0 GeV ^14j . 

We have also checked how the results depend on the intrinsic momentum 
k± by varying its value with 5k±j^ = ±0.1 GeV around our fixed value 
k± = 0.4, which is in the range permitted by the Fermi motion. We have 
repeated the fitting procedure and observed that for both intrinsic momenta, 
k±j^ = k± + Sk±_^ , we are able to find a pair of parameters for which our 
model gives equally good description of data sets as for the central value of 
k± = 0.4. Moreover, we have observed that the value of as(0) parameter for 
all studied intrinsic momenta remains the same but the p±g value is shifted 
for a bigger intrinsic momentum to a higher scale and for a smaller one to a 
lower scale. Therefore, by changing the intrinsic momentum from 0.4 to 0.5 
GeV we can obtain exactly the same best model's parameters set as in |14j 
and the same shape of as presented in Figure 4 from [H]. In that case the 
comparison of the shape of ag in the non-perturbative region of the parton- 
level study are in good agreement with other approaches to modelling non- 
perturbative corrections to inclusive observables with a modified coupling 
in the soft region [T71 [18] . 

The last remark is that using our model as the only non-perturbative 
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ingredient in the simulation, i.e. removing the non-perturbative constituent 
parton masses that usuaUy cut off the parton shower in Herwig++, gives 
a somewhat better description of the data. This lays open the speculation 
that perhaps, in some way, the two approaches could be combined. One 
could for example use our model for the initial-state radiation, and the 
usual model, tuned to describe the final states of e~^e~ annihilation, for 
final-state radiation. 

4. Predictions for LHC and comparison with other approaches 

At the end of this note we would like to compare the results for a trans- 
verse momentum distribution of the Z boson at the LHC energies using the 
nonperturbative gluon emission model and two other approaches: ResBos 
|19| and the Gaussian intrinsic k± extrapolation. But first let us compare 
our prediction of the parton level, marked as the filled histogram in Fig. [3l 
and of the hadron level, dot-dashed blue line. Both histograms, as expected, 
give a consistent extrapolation. 



0.1 I 




5 10 15 20 



p±/GeV 

Fig. 3. Vector boson p± distribution at the LHC. Our model is compared to the 
extrapolation of Gaussian intrinsic fcj^ to LHC energies and the result from ResBos. 



The result from ResBos in Fig. [3] (solid, black) shows a slightly different 
behaviour from our predictions. We predict a slightly more prominent peak 
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and a stronger suppression towards larger transverse momenta. Both com- 
putations match the data well at large transverse momenta as they rely on 
the same hard matrix element contribution for a single hard gluon emission. 
Let us stress the remarkable feature that we predict the same peak position 
with these models which is very important from the experimental point of 
view. This feature is quite understandable as both models are built on the 
same footing: extra emissions of soft gluons. A comparison of ResBos to 
data from experiments at various energies including the experiments E605 
and R209 was done in |20] . 

Furthermore, we see the Herwig++ result from only using intrinsic 
{k±) = 5.7 GeV (dashed, red) as recommended in [12^. This large value 
stems from an extrapolation from lower energy data with the assumption 
that the average k± will depend linearly on ln(M/y^). The peak is seen 
to lie at a considerably higher value of the transverse momentum. It would 
clearly be of interest to have experimental data to distinguish these two 
models of non-perturbative transverse momentum. 

5. Conclusion 

We consider the model based on soft-gluon radiation, much like the 
resummation program ResBos, to have a more meaningful physics input 
than simply extrapolating the Gaussian smearing of a primordial transverse 
momentum. The model implemented in the improved parton shower of 
Herwig++ (release 2.3.1) gives a good (and very similar to the older version 
of Herwig+-|- ) description of data. On the other hand, the fitting procedure 
shows that the best values for the model's parameters are slightly different 
than the previous ones. 
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